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ABSTRACT 

Simnltaneous observations at mnltiple freqnency bands have the potential to 
overcome the fnndamental limitation imposed by the atmospheric propagation 
in mm-VLBI observations. The propagation effects place a severe limit in the 
sensitivity achievable in mm-VLBI, redncing the time over which the signals 
can be coherently combined, and preventing the nse of phase referencing and 
astrometric measnrements. We carried ont simnltaneons observations at 22, 43, 

87 and 130 GHz of a gronp of five AGNs, the weakest of which is ~200 mJy at 
130 GHz, with angular separations ranging from 3.6 to 11 degrees, using the KVN. 

We analysed this data using the Frequency Phase Transfer (FPT) and the Source 
Frequency Phase Referencing (SFPR) techniques, which use the observations at 
a lower frequency to correct those at a higher frequency. The results of the 
analysis provide an empirical demonstration of the increase in the coherence 
times at 130 GHz from a few tens of seconds to about twenty minutes, with FPT, 
and up to many hours with SFPR. Moreover the astrometric analysis provides 
high precision relative position measurements between two frequencies, including, 
for the first time, astrometry at 130 GHz. Finally we demonstrate a method for 
the generalised decomposition of the relative position measurements into absolute 
position shifts for bona fide astrometric registration of the maps of the individual 
sources at multiple frequencies, up to 130 GHz. 

Subject headings: techniques: interferometric, astrometry, radio continuum: galax¬ 
ies (individual: 1803+784, 1807+698, 1842+681, 1928+738 and 2007+777) 


1. Introduction 


Astronomical studies by means of Very Long Baseline Interferometry (VLBI) observa¬ 
tions at cm wavelengths is a well established held, with advanced technological developments 
and analysis tech niques that re sult in superb quality images, including those of very weak 
//Jy sources (e.g. GarrettI 12005 1) and with micro-arcsecond (fias) astrometry measurements 
(e.g. iReid fc Honmal (120141) 1. using phase referencing techniques. This is applied to a wide 
variety of targets and helds of study. 


VLBI at (sub)mm wavelengths {hereafter mm-VLBI) can result in the highest angular 
resolutions achieved in astronomy and has a unique access to emission regions that are in¬ 
accessible with any other approach or at longer wavelengths, because the compact areas of 
interest are often self-absorbed. Therefore it holds the potential to increase our understand¬ 
ing of the physical processes in e.g. Active Galactic Nuclei (AGN) and in the vicinity of 
supermassive Black Holes, and for studies of molecular transitions at high frequencies. 


Nevertheless the applications of mm-VLBI are much less widespread. The observations 
become progressively more challenging as the wavelength gets shorter because of the: lim¬ 
ited telescope surface accuracy and aperture efficiency, receiver system temperatures and 
sensitivity, shorter atmospheric coherence times and that sources are intrinsically weaker in 
general. Moreover phase referencing techniques,which are routinely used in cm-VLBI, fail to 
work bey ond 43 GHz (excluding a single case at 86 GHz with the Very Long Baseline Array 
(VLBA) ( Boreas fc Rioia 20021) 1. 


Gontinuous development and technical improvements have led to a sustai ned increase of 
the hig h frequency threshold for VLBI observations in the last two decades fsee lKrichbanm et al. 
(120141) for a review). Regular observations up to 86 GHz are being carried out with well es¬ 
tablished networks such as the VLBA and the Global Millimeter VLBI Array (GMVA), more 
recently up to 130 GHz with the Korea n VLBI Network (KVN ) and ad-hoc observations at 
the highest frequencies up to 240 GHz (iDoeleman et al.l 120081) . The held of mm-VLBI will 
greatly beneht from the a rrival of the Phased-up Atacama Large Millimeter Array (ALMA) 
(iMatthews fc Grew! 120151) for joint VLBI observations. 


In this paper we will focus on two aspects that limit the potential of mm-VLBI observa¬ 
tions: 1) achieving improved sensitivity through increased coherence times, to increase the 
number of targets; 2) achieving astrometry, and in particular for “bona hde” registration of 
images at multiple frequency bands, to reveal the physical processes in a number of helds of 
astronomy. 


For example, for AGN studies, maps of the Spectral Index or Rotation Measure across 
the source, at mm-wavelengths, provide crucial insights into the development of the magnetic 
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field strength and particle densities as the jet exits the core region and extends down the 
out-flow. However, astrometric map registration is crucial to make a reliable measurement 
and to form meaningful interpretation. There are a number of methods in AGN studies 


which can be used to align images at multiple frequencies, as discussed in iHovatta et ah 


(I 2 OI 4 J) . They argued that any results derived without accurate astrometric registration are 
questionable in the vicinity of the core, which is the most interesting region in mm-VLBI. 
Also, in studies of the maser emission from the molecular species that exist in circumstellar 
envelopes (CSEs) and star forming regions (SFR), the comparison of the locations of the 
different species of maser emission can be inverted to reveal the physical conditions as a 
function of the dista nce from the central star pumping the masers (e.g. see discussion in 
Reid &: Moran f IQSsh i. In both cases this would allow one to fully understand the flow of 
material and energy in stellar environments during the formation and the evolution of stars. 
Traditionally there has been no other mechanism other than phase referencing to accurately 
astrometrically register the maps at the different bands. 


The Korean VLBI Network (KVN) (iKim et al.l |2004J: Ihee et al.ll201# is the hrst ded¬ 
icated mm-VLBI array and addresses one of the fundamental limitations of the held, the 
atmospheric stability. It currently consists of three antennas operated by Korea Astronomy 
and Space Science Institute (KASI), spread across South Korea, located in the campus of the 
Universities of Yonsei and Ulsan in main land and on Jeju island. The observing frequencies 
are centred at 22, 43, 87 and 130 GHz. The baseline lengths between the antennas range 
between 300 and 500 km, which provide a spatial resolution ^1 m as at the highest frequency 


band. The innovative multi-band receiver flHan et al.ll2008l. l2013l) of KVN is designed to mit¬ 
igate the atmospheric propagation effects using simultaneous observations at multiple bands. 
The KVN combined with the Frequency Phas e Transfer (FPT) and Source Frequency Phase 


Referencing (SFPR) da ta analysis techniques (jPodson fc Rioiall2009l: iRioia fc Dodsonll2011 


Rioia et al.ll201ll. I2014J ) (see also references therein) allows an effective increase of the coher¬ 
ence time, well beyond that imposed by tropospheric fluctuations, as well as high precision 
astrometric measurements, respectively, even at the highest frequencies. We know of no 
demonstrated upper frequency limit and the methods would be expected to work as long as 
the tropospheric propagation effects were non-dispersive. 


In this context, successful tests are on-going with ALMA flFomalont et al.l 1201411 at fre- 
quencies as high a s 650 GHz (where this is known as the Band-to-Band (B2B) mode). In 
Rioia et al.l (1201411 we presented results of SFPR astrometric measurements with KVN at 
22 and 43 GHz for continuum sources, along with a detailed cornparati ve study using fast 
frequency switching observations with the VLBA. iDodson et al.l (1201411 presented the ap¬ 
plication of SFPR to spectral line studies for astrometric registration of the H 2 O and SiO 
maser maps, at 22 and 43 GHz, respectively, in GSEs. In this paper we extend the astromet- 












































4 


ric measurements to all the four frequency bands supported by the KVN, up to 130 GHz, 
and quantify the increase in coherence time. The paper layout is as follows: The simulta¬ 
neous multi-frequency observations at four-bands are presented in Section 2; a description 
of the analysis carried out to obtain the maps, astrometric measurements and “bona hde” 
astrometric registration of multi-frequency images is in Section 3; the results are presented 
in Section 4 and a discussion of the results in Section 5. 


2. Observations 


In March 5, 2014, we carried out simultaneous observations at four frequencies, i.e. 22, 
43, 87 and 130GHz (also known as K, Q, W and D bands, respectively), with the three 
antennas of the KVN, towards 5 AGN target sources, for a total duration of 9 hours. 

The recording consisted of 4 consecutive 16 MHz intermediate frequency (IF) sub-bands 
at each frequency. The lower edges of the first IFs at Q, W and D bands were selected to be 
multiples of that at K band, those being 21.65, 43.30, 86.60 and 129.90 GHz. Having integer 
frequency ratios is important for the successful application of tropospheric compensation 
techniques using multi-frequ ency observations. The correlation was done with the DiFX 
correlator flDeller et al.ll201ll) with 1 second averaging and a spectral resolution of 64 channels 
per IF. 


The target sources were selected from the 86-GHz VLBI catalog (ILee et al.ll2008[) based 
on two criteria: to have strong detections at W band, and angular separations in the sky 
ranging from a few to many degrees. Figure [T] shows the distribution in the sky of the 
five selected sources (1803-1-784, 1807-1-698, 1842-1-681, 1928-1-738 and 2007-1-777) along with 
their angular separations, which range between ~ 3.6° and 11°. None of the sources had 
been observed previously with VLBI at 130 GHz. The observations consisted of ~3 minutes 
long scans alternating between the sources in each of the two triangles shown in Figure 
m and between the triangles ca. every hour and a half, at the four bands simultaneously. 
Alternating between multiple sources allowed us to develop a strategy to decompose the 
relative astrometric measurements into single source position shifts that allow, for example, 
the registration of the images at different frequencies. 


3. Methods 


In this section we describe the various mapping and astrometric analyses carried out: 








5 



Fig. 1.— Sky distribution of the five AGN sources in this study. The source pairs in the 
SFPR analysis are connected with a line, and have angular separations ranging from 3.6° 
to 11°. The triangles connect sources that were observed in a ~1.5 hour block, alternating 
between the two triangles. 
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3.1. Hybrid Maps at the Four KVN Frequency Bands. 


19951). 


_We followed standard procedures for imaging VLBI datasets using AIPS flDiamond 
One of the major challenges of imaging KVN observations arises from the small 


number of antennas, which prevents the application of amplitude self-calibration techniques 
to derive amplitude gain corrections, as this requires a minimum of four antennas. 

At the lower frequencies (K and Q bands), the system temperature measurements, along 
with regul ar sky-dipping, have been shown to provide a good estimate of the system per¬ 
formance (ILee et al.ll201411 . However, at the higher frequencies (W and D bands) significant 
discrepancies can be expected. Hence, we have attempted to use the observations of the 
hve target sources to estimate global (i.e. for all sources) amplitude gain correction factors 
that should rescale the nominal calibration information, at each band. We assumed a point 
source model of arbitrary flux (for all sources) and calculated the normalized amplitude gain 
corrections, for each individual source and at each frequency band. The individual gain 
estimates for all sources at a given band showed a good agreement, as shown in Figure [21 
this supports the validity of the assumption of point source structure with KVN resolutions 
at all bands. At each band, the values for all sources were merged and smoothed together, 
thereby further suppressing any individual source structure contributions to the estimated 
gain amplitudes, except for one of the weakest sources, 2007-1-777, at the highest frequency, 
130GHz, which has noisy solutions. The resultant amplitude gain correction factors were 
applied to the corresp onding datasets using AIPS. The hybrid maps, made with difmap 
([Shepherd et al.lll994l) . are presented in Section [H 


3.2. SFPR Maps and Astrometric Analysis 


We carried out SFPR astrometric analysis of the four-band multi- frequency KVN dataset 


using AIPS. The details of SFPR analysis are presented elsewher e (iDodson fc Rioial 12009 


Rioia fc DodsorJ[20Il[: [Rioia et al.[[20Il[. [20I4J: [Dodson et al.[[20I4j) and basically consists of 
two calibration steps. In a hrst step, the observations at the higher frequency bands (z/^‘®^) 
are calibrated using the simultaneous observations at a lower frequency band (p*°"), for each 
source. This is done for all frequency pairs which have an integer frequency ratio R (with 
R = yis^), which in turn is used to scale up the phase-calibration solutions at the lower 
frequency. This dual frequency calibration step eliminates the common non-dispersive resid¬ 
ual errors (e.g. tropospheric propagation effects and inaccurate coordinates) in the complex 
visibilitity output of the correlator, providing an increased signal coherence at the higher 
frequency. We have dubbed this step “Frequency Phase Transfer” (FPT) and the outcome 
FPT-ed visibilities. The second step of the calibration removes the remaining dispersive 
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Fig. 2.— Superimposed normalized amplitude gain correction factors to the nominal KVN 
calibration, versus time, derived assuming point source models for each source (small sym¬ 
bols: ‘x’, circles, squares, triangles and diamonds for 1803-1-784, 1807-1-698, 1842-1-681, 
1928-1-738 and 2007-1-777 respectively), at all bands: a) 22 GHz, b) 43 GHz, c) 87 GHz and 
d) 130 GHz. At each band, the corrections from the individual sources were merged and 
smoothed to provide the final set of amplitude corrections, which were globally applied (large 
black crosses). At the lower frequencies (i.e. 22 GHz and 43 GHz) the estimated corrections 
were close to unity. On the other hand, at the higher frequencies the corrections spanned 
a greater range both above and below the nominal value, with long periods where the gain 
corrections are less than unity. That implies significantly better than 100% efficiency, which 
is unlikely, therefore these gains were adjusted so that the maximum efficiency correction was 
100%. These occurred at the highest elevations where one would expect the least deviation 
from the nominal gain performance. 




residual errors (i.e. instrumental and ionospheric propagation effects) using the interleaving 
observations of another source. This two-step calibration retains the astrometric signature 
of any source position shifts between the two frequencies in the interferometric phase. The 
resultant calibrated visibilities, for a given frequency pair and a source pair, are dubbed 
SFPR-ed visibilities. The Fourier transformation of the SFPR dataset is the SFPR map, 
which conveys a bona fide astrometric measurement of the relative separation or shift be¬ 
tween the position of the reference points in the images at the two frequencies, for the two 
sources. If the reference points are selected to be the “core” components this shift corre¬ 
sponds to the “core-shift” phenomenae, or change o f position at the bas e of th e jet, due 
to fre quency dependent opacity effects as defined by iBlandford fc Konigll fjl979[) : iLobanov 
fll998l) . for AGN observations. In general, the KVN resolution will result in structure blend¬ 
ing of the “core” and other jet components into a single unresolved component in the maps, 
therefore this shift gives the relative angular separation between the centroids of the bright¬ 
ness distributions at the two frequencies, for the two sources. Through-out this text we are 
using the term “position-shift” (rather than “core-shift”) to include all cases. For a detailed 
description of structure blending effects in astrometric analysis see Rioja et ah 2014. 


In this paper, the SFPR analysis comprises of five pairs of frequency bands, with values 
of R ranging from 2 to 6, and six pairs of sources, with angular separations between 3.6° 
and 11°. The frequency band pairs along with the corresponding scaling factors for the 
lower frequency (R) in parenthesis are: K—)-Q (x2), K—)-W (x4), K—)-D(x6), Q—)-W(x2) 
and Q—)-D(x3). The six pairs of sources are shown in Figured] with connecting lines. 


3.3. Astrometric Registration of Images Across Frequency Bands 

The outcome of the SFPR analysis are relative astrometric measurements and, just 
like the measurements from phase referencing analysis, there is an inherent ambiguity as 
to what are the individual contributions arising from each of the two sources; that is, the 
solutions are degenerate. Increasing the number of sources sets stronger constraints in the 
disentangling into individual contributions, as the additional constraints reduces, but does 
not break the degeneracy. Our observations comprise of multiple pairs of sources for this 
purpose. Starting with the pairwise astrometric SFPR measurements, we estimated the 
single source position-shifts, using the Singular Value Decomposition (SVD) as the linear 
least squares minimization method, for each frequency pair. 

No matter how many combinations of sources one measures there will always remain 
an ambiguity of the global (i.e. common for all sources) absolute correction. That is, for 
example, in the case that all the sources have identical position-shifts it would leave no 
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signature in the measurements. Therefore, once we have decomposed the measurements 
into the contributions from the individual sources, we still need to find the global absolute 
correction. We have additional information which allows us to do this in most cases, namely, 
that a frequency dependent source position shift is expected to be aligned with the direction 
of the jet in the map of the source. This expectation applies for both types of position 
shifts described above, i.e. core-shifts arising from opacity effects and/or centroid shifts 
arising from structure blending, regardless of its nature. This approach will fail for the cases 
when the position-shifts and the jet axis do not coincide, which would be unexpected, or 
when all the sources have similar jet directions and therefore there is no clear best global 
correction to be determined. The group of sources in our ob servations shows a wide range of 
jet directions in the high resolution VLBI MOJAVE maps flLister fc Homan!l2005h . Hence, 
by adding the constraint that the position-shift direction must align with the up-stream 
jet direction we can determine both the appropriate global correction and unambiguous 
individual source position-shifts (also called absolute position shifts hereafter). The latter 
are also the shifts required for a bona-fide astrometric registration of the images at the four 
observed frequencies, for all the sources. 


4. Results 


4.1. Hybrid Maps at the Four KVN Frequency Bands. 


Figure [3] shows self-calibration images for the hve AGN sources in the multi-frequency 
KVN observations, including the hrst images of these sources at 130 GHz. The visibility 
datasets were modelfitted and imaged and show little divergence from point-source core 
dominated images; in some cases there appears to be some elongation aligned wit h the jet 
direction as seen in higher resolution maps from MOJAVE (iLister fc Homan! 120051) . 


Table d] lists the total flux values as measured from the maps, at all frequency bands. 
It should be noted that the absolute flux values might suffer from the lack of absolute cali¬ 
bration, especially at the highest frequency band, 130 GHz (see section 3.1 for a description 
of the amplitude calibration). 


It is worth emphasizing that not all sources had direct detections, i.e. within the atmo¬ 
spheric coherence time, and for those that didn’t, we benefited from the extended coherence 
time resulting from a previous trans-frequency FPT analysis calibration. Note that remain¬ 
ing dispersive residuals prevent making a map after solely FPT calibration; nevertheless the 
FPT analysis conditions the dataset at and allows for a self-calibration analysis us¬ 
ing much longer phase solution time intervals, hence enabling the detection of sources that 
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would not be detected otherwise (i.e. within the atmospheric coherence time interval). The 
resulting maps are therefore self-calibration maps, and have no astrometry information. This 
procedure enabled the imaging of 2007-1-777 and 1842-1-681 at 130 GHz, which were too weak 
for direct detections. 



22 GHz 

43 GHz 

87 GHz 

130 GHz 

1803+784 

1.7 

1.2 

1.0 

0.6 

1807+698 

1.1 

0.9 

0.9 

0.5 

1842+681 

0.5 

0.4 

0.4 

0.2 

1928+738 

3.7 

3.2 

3.3 

1.6 

2007+777 

0.6 

0.4 

0.4 

0.3 


Table 1: Total Source flux, in Jy, for the five sources and at the four frequency bands, 
measured from the self-calibrated maps in Figure |3l 


4.2. Increased Coherence Time for mm-VLBI 

The rapid changes in the observed interferometric phases introduced by the tropospheric 
propagation effects set a severe limit on the coherence time for integration of the signal in 
observations at high frequencies, and therefore the sensitivity of those observations. A direct 
consequence of the effective tropospheric compensation achieved from simultaneous dual 
frequency observations is an increased coherence time and therefore sensitivity. This can be 
visually appreciated in the FPT and SFPR-ed phases shown in Figures 0] and 0 respectively. 
Figure m shows the FPT-ed calibrated visibility phases at for the hve frequency pairs 
(j^iow j/high^ with R integer, shown in separate plots, in our observations. Note that in all 
cases, the calibration applied has been derived from a different frequency band and scaled 
with the corresponding factor R. In all cases, the compensation of the fast tropospheric 
fluctuations results in a much higher degree of coherence, compared to the raw output of the 
correlator. Figure 0] shows the SFPR-ed visibility phases for a subset of frequency pairs and 
source pairs (using 1803-1-784 as reference) in these observations, which are representative 
of the final products of the SFPR analysis. It is immediately obvious that the remaining 
dispersive residual phase variation in the FPT-ed visibilities has been compensated for in 
the SFPR visibilities. 

In order to quantify the increase in the coherence time, we have carried out a comparative 
study of the coherence times achieved with FPT and SFPR at 130 GHz. This is the highest 
frequency in our observations, where the propagation effects are most severe. To perform 
these tests we use the AIPS task CALIB on FPT-ed and SFPR-ed calibrated datasets of 
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22 GHz 43 GHz 87 GHz 130 GHz 



10 0 -10 10 0 -10 5 0 -5 5 0 -5 pq 



10 0 -10 10 0 -10 5 0 -5 5 0 -5 



10 0 -10 10 0 -10 5 0 -5 5 0 -5 

RA offset (mas) 


Fig. 3.— Self-calibrated images of the five AGNs at the four KVN frequency bands. The 
contours in all cases start at ±50mJy and double for each subsequent contour line. The 
image sizes are ±10 mas at K and Q bands, and ±5 mas at W and D bands. The beam size 
is indicated by the elipse at the bottom left of the image. 
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Time (hours) Time (hours) 

Fig. 4.— Superimposed FPT-ed residual visibility phases for all the sources in this study, 
in separate plots for each of the hve frequency pairs —)• that have an integer 

frequency ratio R. Note that the FPTed visibility phases are the correlator output at 
solely calibrated with the scaled up phase solutions at Upper row: Corresponds to 

frequency pairs with = 22 GHz (K band). From left to right, K—)-Q (R=2), K—)-W 
(R=4) and K—)-D (R=6). Lower row: Corresponds to frequency pairs with = 43 GHz 
(Q band). From left to right, Q—j-W (R=2) and Q—)-D (R=3). Different colours and symbols 
correspond to observations of different sources: red (‘x’) for 1803+784, green (‘circle’) for 
1807+698, orange (‘+’) for 1842+681, dark blue (‘square’) for 1928+738 and cyan (‘triangle’) 
for 2007+777. The enhanced coherence time is immediately obvious as is the agreement 
between the phases for the different sources, for each frequency pair. The rightmost plots 
underline the benehts from using the scaled calibration values derived from = 43 GHz 
calibration for the 130-GHz data, as compared to = 22 GHz. 
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a) 1928+738: K-W 


b) 2007+777: Q-W c) 1807+698: Q-D 


d) 1842+681: Q-D 
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Fig. 5.— Outcomes of SFPR astrometric analysis: Upper Row: SFPR-ed residual visibility 
phases for a subset of source and frequency pairs from the comprehensive analysis presented 
in this paper. The target source and target frequency band (i.e. along with the 

angular separation to the reference source (in all cases 1803+784), the reference frequency 
band (i.e. and the frequency ratio R, in parenthesis, are specihed for each plot. From 

left to right: 1928+738 at W-band (6.8° apart on the sky, K-band, R=4); 2007+777 at W- 
band (6.3° apart, Q-band, R=2); 1807+698 at D-band (8.6° apart, Q-band, R=3); 1842+681 
at D-band (10.7° apart, Q-band, R=3). Lower Row: SFPR-ed maps resulting from Fourier 
transformation of the SFPR-ed visibility phases directly above. From left to right: 1928+738 
at 87GHz (W-band), 2007+777 at 87 GHz; 1807+698 and 1842+681 at 130 GHz (D-band). 
Peak fluxes are 2 Jy beam“^, 266 mJy beam“^, 415 mJy beam“^ and 216 mJy beam“^, 
respectively. The contour levels in the maps start from 0.75% of the corresponding peak 
fluxes, respectively, and doubling thereafter in all cases. Each map includes a negative 
contour level at the same percent level of the peak flux as the hrst positive one. The beam 
size is indicated at the bottom left of the image. 
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Fig. 6.— Plot of the Fractional Flux Recovery quantity (see text for details) versus length of 
(self-calibration) phase solution time interval, in minutes, for comparative coherence studies 
at 130 GHz between FPT (blue dashed line, with box symbols) and SFPR (red solid line, 
with circle symbols) calibration. The coherence time corresponds to a fractional flux recovery 
of ~ 0.6. For this study we used KVN observations of 1842-1-681 at = 130 GHz (D- 
band) calibrated with scaled up phase solutions at = 43 GHz (Q-band) for FPT, plus 
1803-1-784 as the second source for SFPR analysis. 
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1842+681 at = 130 GHz, with = 43 GHz and 1803+784 as the reference sonrce, 
with a series of phase solntion time intervals ranging from 0.5 to 480 minntes. In each case 
the phase solntions are applied and the calibrated visibility data Fourier inverted to produce 
a map. We use the fractional peak recovered flux quantity, defined as the ratio of the peak 
flux in this map and that from self-calibrated maps, as a measure of remaining phase errors 
in the analysis. Figure |6] shows that the fractional peak recovered flux values in the maps 
decrease with increasing temporal solution intervals, as expected. The coherence time is 
defined as the solution interval at which the peak flux recovery is 60%, being equivalent to 
the rms residual phase being equal to 1 radian. Our analysis show that the coherence time at 
130-GHz is ~20 minutes with FPT calibration. With SFPR calibration there is practically 
no limit in the coherence time; we could integrate up to 8-hours, the whole duration of the 
experiment, with a mere 20% loss of peak flux. Note that the tropospheric coherence time 
at 130-GHz is some tens of seconds and that neither 1842+681 nor 2007+777 have direct 
detections at this frequency. 


4.3. SFPR Maps and Astrometry at 22, 43, 87 and 130-GHz. 


The hnal outcome of the SFPR analysis is a SFPR map which conveys the astrometric 
information. Figure E] shows a subset of the SFPR-ed maps obtained in the comprehensive 
analysis of KVN observations; they are the fourier transform of the SFPR visibilities directly 
above in the same figure. The offset of the peak of brightness with respect to the center of 
the maps is a measurement of the relative position shift between the two frequency bands, 
for the two sources. The complete results from the SFPR astrometric analysis, comprising 
of the hve frequency pairs and six source pairs are summarized in Table [2J Table [2] lists the 
right ascension and declination offsets of the peak of brightness from the center of the SFPR 
maps, as measured with AIPS task MAXFIT; the dynamic range of the map; the estimates 
of rms SFPR phase e rrors arising from the different contributions following formulae in 
Rioia fc DodsoiJ (1201 ll) . along with their quadratic sum; the last columns are the estimated 
astrometric errors as described in Section 4.5. Note that the dual frequency calibration 
provides a perfect compensation for the non-dispersive phase errors (i.e. 0geo*^, and ^^tfp > 
which stand for geometric, dynamic and static tropospheric errors, respectively), that the 
thermal noise term (0]))^) and the dynamic ionospheric errors are rarely significant 

and that the phase errors are dominated by the contribution from the static ionosphere 
contribution latter is largest for the frequency pairs with = 22 GHz and 

increases with the angular separation between the sources. 
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‘r<PS 

(°) 

n 

^Aa cos 5 

(/ ias ) 

O'AS 

(/ ias ) 




1803+784 / 1807+698 ( 8 . 65 '’ 

apart) 






K^Q 

86 

7 

390 

1.6 

0.(^6 

0.0 

0.0 

4.1 

24.6 

25.0 

88 

127 

K^W 

81 

76 

320 

1.9 

0.0 

0.0 

0.0 

10.3 

61.5 

62.4 

110 

159 

K^D 

129 

132 

23 

27.1 

0.0 

0.0 

0.0 

16.1 

95.7 

100.7 

118 

171 

Q^W 

-11 

65 

360 

1.7 

0.0 

0.0 

0.0 

2.1 

12.3 

12.6 

22 

32 

Q^D 

10 

82 

138 

4.5 

0.0 

0.0 

0.0 

3.7 

21.9 

22.6 

27 

38 




1803+784 / 1842+681 ( 10 . 70 ' 

’ apart) 






K^Q 

-200 

75 

360 

1.7 

0.0 

0.0 

0.0 

4.4 

30.4 

30.8 

109 

157 

K^W 

-364 

169 

200 

3.1 

0.0 

0.0 

0.0 

11.1 

76.1 

76.9 

136 

196 

K^D 

-379 

241 

40 

15.6 

0.0 

0.0 

0.0 

17.2 

118.3 

120.6 

142 

205 

Q^W 

-192 

110 

240 

2.6 

0.0 

0.0 

0.0 

2.2 

15.2 

15.6 

27 

40 

Q^D 

-235 

130 

57 

10.9 

0.0 

0.0 

0.0 

3.9 

27.0 

29.4 

35 

50 




1803+784 / 1928+738 ( 6 . 79 '’ 

apart) 






K^Q 

-86 

70 

830 

0.8 

0.0 

0.0 

0.0 

3.9 

19.3 

19.7 

69 

100 

K^W 

-241 

150 

380 

1.6 

0.0 

0.0 

0.0 

9.7 

48.3 

49.3 

87 

125 

K^D 

-307 

210 

100 

6.2 

0.0 

0.0 

0.0 

15.0 

75.1 

76.8 

90 

130 

Q^W 

-156 

64 

960 

0.6 

0.0 

0.0 

0.0 

1.9 

9.7 

9.9 

17 

25 

Q^D 

-165 

84 

150 

4.2 

0.0 

0.0 

0.0 

3.4 

17.2 

18.0 

21 

31 




1803+784 / 2007+777 ( 6 . 34 '’ 

apart) 






K^Q 

-43 

43 

390 

1.6 

0.0 

0.0 

0.0 

3.8 

18.0 

18.5 

65 

94 

K^W 

0 

18 

70 

8.9 

0.0 

0.0 

0.0 

9.5 

45.1 

46.9 

83 

119 

K^D 

-4 

45 

69 

9.0 

0.0 

0.0 

0.0 

14.8 

70.1 

72.2 

85 

123 

Q^W 

-41 

15 

540 

1.2 

0.0 

0.0 

0.0 

1.9 

9.0 

9.3 

16 

24 

Q^D 

-42 

50 

85 

7.3 

0.0 

0.0 

0.0 

3.4 

16.0 

17.9 

21 

30 




1807+698 / 1842+681 ( 3 . 60 '’ 

apart) 






K^Q 

-276 

73 

376 

1.7 

0.0 

0.0 

0.0 

3.4 

10.2 

10.9 

38 

56 

K^W 

-476 

123 

227 

2.7 

0.0 

0.0 

0.0 

8.5 

25.6 

27.1 

48 

69 

K^D 

-549 

156 

67 

9.3 

0.0 

0.0 

0.0 

13.2 

39.8 

43.0 

50 

73 

Q^W 

-185 

47 

133 

4.7 

0.0 

0.0 

0.0 

1.7 

5.1 

7.1 

13 

18 

Q^D 

-230 

45 

143 

4.4 

0.0 

0.0 

0.0 

3.0 

9.1 

10.5 

12 

18 




1928+738 / 2007+777 ( 4 . 52 '’ 

apart) 






K^Q 

125 

-92 

193 

1.9 

0.0 

0.0 

0.0 

3.5 

12.9 

13.5 

47 

69 

K^W 

181 

-205 

162 

2.2 

0.0 

0.0 

0.0 

8.8 

32.1 

33.4 

59 

85 

K^D 

192 

-257 

85 

4.2 

0.0 

0.0 

0.0 

13.7 

50.0 

52.0 

61 

88 

Q^W 

69 

-76 

73 

4.9 

0.0 

0.0 

0.0 

1.8 

6.4 

8.3 

15 

21 

Q^D 

210 

-105 

70 

5.1 

0.0 

0.0 

0.0 

3.1 

11.4 

12.9 

15 

22 


Table 2: Summary of the measurements from the SFPR astrometric analysis presented in this paper, along 
with the error estimates, for the five frequency pairs (column 1) and for the six source pairs (separated by 
horizontal lines). The relative astrometric offsets and the dynamic ranges (columns 2-3 and 4, respectively) 
are measured from the SFPR maps. A list of the estimated error contributions, per baseline, is provided: the 
thermal errors (column 5) are estimated from the dynamic range (see text); the geometric and propagation 
media contribution errors (columns 6-10) are estimated using the formulae in iRioia fc DodsonI (|201111 . for 
typical parameter uncertainties of the tropospheric zenith path delay and the TEC equal to 3 cm and 3 
TECU, respectively, source angular separations as listed, simultaneous multi-band observing (Tj)^^ = 0) and 
source switching cycle of Tswt = 450 seconds. Column 11 is the quadratic sum of the forementioned errors. 
Columns 12 and 13 are the errors of the SFPR astrometric measurements, in right ascension and declination. 
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Source Name 

1803+784t 

1807+698t 

1842+681 

1928+738^ 

2007+777 

Jet PA (deg) 

91.5T3.5 

97.5T2.5 

-135+5 

-160+10 

90+5 


able 3: Jet position angle (PA) and error range from high resolution maos in iHovatta et ah 


(120141) where marked with P otherwise measured from the MOJAVE maps flLister &: Homan 


200511 with an estimated error of 5°. 


Freq. Pairs 

K^Q 

K^W 

K^D 

Q^W 

Q^D 

Global Correction (pas) 

-42,-36 

-82,-76 

-88,-100 

-48,-326 

-82,-34 

Errors (pas) 

+14,+4 

+29,+10 

+39,+13 

+10,+6 

+14,+9 


Table 4: List of the global astrometric correction vectors, in pas on the sky of right ascension 
and declination, that result in the best alignment of the jet directions (listed in Table [3]) and 
the individual source frequency dependent position shifts, for the five frequency pairs. The 
errors in these corrections are given below. 
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4.4. Decomposition of relative astrometric measurements into individual 

source position-shifts 


Figure [7] shows polar plots of the pairwise astrometric measurements listed in Table 
[21 for the four source pairs involving 1803+784, and for the hve frequency pairs. These 
are direct outcomes of the SFPR analysis. For each source pair, the measurements are the 
combined position shift contributions from both sources between the two frequencies, for each 
frequency pair and therefore are expected to show little directional coherence, except when 
one source has a dominant position-shift (e.g. the plots involving 1842+681 or 1928+738). 


We have used SVD to decompose the pairwise position shifts into single source frequency 
dependent position shifts, albeit with degeneracies included. Those are shown in Figure |8l 
where one can appreciate an improved agreement between the directions of the position- 
shifts for each source, although those are not well aligned with the jet direction in high 
resolution maps. The jet directions in the high resolution maps for the hve AGN sources 
flHovatta et ah!I2014J: [Lister fc Homanll2005[) are also shown in Figure [Hfe and their values 
are listed in Table [3l We used the expectation of alignment to break the degeneracy, by 
hnding the best global correction (i.e. common for all sources) through a grid search of a 
few hundred micro-arcseconds around the SVD solutions, for each frequency pair. Figure [9] 
shows the degree of alignment as a function of grid position for the K—)-Q dataset. Table [H 
lists the global corrections that were found to best align the SVD minimized single source 
position shifts to the jet directions, for the hve frequency pairs. Finally, Table [5] lists the 
resultant absolute single source position shifts corresponding to the hve frequency pairs, for 
the hve sources. Figure [10] shows polar plots of these single source frequency dependent 
absolute position shifts listed in Table El which display a tight agreement between the hve 
frequency pairs and are well aligned with the jet directions for each source. 
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3)1803+784/1807+698 

0 



b)1803+784/2007+777 
0 



c)1803+784/1928+738 d)1803+784/1842+681 

0 0 




Fig. 7.— Polar plots of the SFPR astrometric measurements in Table [21 corresponding to 
the four pairs of sources with 1803+784 (plots a-d, with 1807+698, 1842+681, 1928+738 and 
2007+777, respectively) at the hve pairs of frequencies (shown in different colors): K—)-Q 
(red), K—)-W (blue), K—)-D (black), Q—)-W (green) and Q—)-D (cyan). The vectors corre¬ 
spond to the relative position shifts for the two sources between two frequencies, for the hve 
frequency pairs. In the polar plots the position angles are shown outside the largest circle 
and are 0 and 90 degrees towards North and East, respectively, and the magnitude units, as 
specihed in the concentric circles, are in ^uas. 
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a) 1803+784 


0 



d) 1928+738 
0 



b) 1807+698 


0 



e) 2007+777 
0 



c) 1842+681 


0 



f) Jet Directions 


0 



Fig. 8.— (a-e) Polar plots of the degenerate decomposed single source position shifts between 
two frequencies, for the hve frequency pairs and the hve AGNs. These are shown in different 
colors (as in Figure [7]): K—)-Q (red), K—)-W (blue), K—)-D (black), Q—)-W (green) and Q—)-D 
(cyan). They have been derived from the SFPR pair-wise measurements using SVD. Each 
vector represents the best linear least squares minimised position shift, and includes a global 
(i.e. common for all sources) degenerate offset for each frequency band. The axes of the 
polar plots are as in Figure [71 (f) Polar plot of the jet directions as appear in high resolution 
maps of the hve AGNs, showing the poor alignment with the position shifts in (a-e). The 
vector colours for the jet directions are: red for 1803-1-784, green for 1807-1-698, orange for 
1842+681, blue for 1928+738 and cyan for 2007+777. 
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Fig. 9.— The combined misalignment for all sources (Z-axis) between the jet directions in 
Table 3 and the frequency dependent position shift vectors of the same source, evaluated on 
a grid of global (i.e. common for all sources) corrections to the degenerate SVD position 
shift vectors shown in Figure 8, for the frequency pair 22/43 GHz (K—)-Q). The grid ranges 
from ±200pas in right ascension and declination (X and Y-axis, respectively). The global 
misalignment is shown in arbitrary units of the absolute sum of the differences between the 
complex direction vectors of jets and modihed position shift measurements for all sources. 
The best global alignment between the directions of the jet and the position shift corresponds 
to a correction of -42 and -36 pas, in right ascension and declination, respectively, to the 
SVD solution. A similar procedure was carried out for each frequency pair; the results of 
the complete analysis are listed in Table 01 
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a) 1803+784 


0 



d) 1928+738 
0 



b) 1807+698 


0 



e) 2007+777 

0 



c) 1842+681 


0 



f) Jet Directions 


0 



Fig. 10.— (a-e) Polar plots of the decomposed absolute single source position shifts between 
two frequencies, for the hve frequency pairs and the hve AGNs. These are shown in different 
colors (as in Figure [7]): K—)-Q (red), K—)-W (blue), K—)-D (black), Q—)-W (green) and Q—)-D 
(cyan). They have been derived from the SFPR pairwise measurements, using SVD plus 
the alignment constraint between the jet and the position shift directions, to break the 
degeneracy. The axes of the polar plots are as in Figure [71 (f) Polar plot of the jet directions 
as for Figure |8l showing the good alignment between the jet directions and that of the 
decomposed frequency dependent position-shifts, for each source. Thus we have obtained 
absolute position shifts for the individual sources from the original pairwise measurements. 




























23 


1 

2 3 4 5 6 

Absolute Position Shift Registration Error 

Freq. 

Pair 

Aa cos5 
(/xas) 

(/ias) 

PA (^AacosS 
(deg) (/xas) 

ctas 

(//as) 

K^Q 

-92 

-4 

1803+784 

92 

31 

41 

K^W 

-200 

0 

89 

43 

51 

K^D 

-225 

19 

84 

50 

54 

Q^W 

-132 

18 

82 

11 

12 

Q^D 

-173 

30 

80 

15 

15 

K^Q 

-164 

-3 

1807+698 

91 

24 

29 

K^W 

-272 

-72 

104 

36 

37 

K^D 

-324 

-93 

106 

43 

39 

Q^W 

-116 

-56 

115 

10 

9 

Q^D 

-167 

-59 

109 

14 

12 

K^Q 

109 

-78 

1842+681 

-125 

66 

93 

K^W 

183 

-182 

-134 

85 

117 

K^D 

188 

-235 

-141 

91 

122 

Q^W 

64 

-97 

-146 

17 

23 

Q^D 

61 

-102 

-148 

24 

30 

K^Q 

17 

-77 

1928+738 

-167 

32 

43 

K^W 

46 

-145 

-162 

45 

54 

K^D 

85 

-204 

-157 

52 

57 

Q^W 

16 

-42 

-158 

12 

12 

Q^D 

27 

-54 

-153 

15 

16 

K^Q 

-78 

-16 

2007+777 

101 

18 

20 

K^W 

-167 

20 

82 

31 

26 

K^D 

-163 

13 

85 

39 

28 

Q^W 

-71 

18 

75 

10 

9 

Q^D 

-157 

15 

84 

13 

10 


Table 5: List of absolute single source position shifts between two frequencies, for the five frequency 
pairs and for the five AGNs in this study, as plotted in Figure [TOl They have been derived from the 
SFPR pairwise measurements, using SVD plus the alignment constraint between the jet and the 
position shift directions, to decompose into single source contributions and to break the degeneracy. 
Column 1 are the frequency pairs; Columns 2 and 3 are the Right Ascension and Declination of 
the position shifts, respectively, with the corresponding position shift direction (PA) in Column 4. 
The errors (Columns 5 and 6) include all random and systematic contributions. The position shift 
direction (PA) should be compared with the jet position angles given in Table El 











24 


4.5. Astrometrical Error Analysis 


We have carried out a comprehensive error analysis to estimate the propagation of 
random and systematic error contributions in the SFPR analysis, along with those from 
the SVD and global-shift minimization analysis, into the frequency dependent position shift 
astrometric measurements for each source. 


For the SFPR error analysis, we have used the formulae in iRioia fc Dodson! (120111) to es¬ 
timate the residual phase errors arising from typical parameter uncertainties in the ‘a priori’ 
models for the propagation medium and the geometry contributions. The estimated values 
per baseline are listed in Table [H for the geometry and for the dynamic and static 

components of both the troposphere {cr(pdtrp cy(f)ll^p, respectively) and the ionosphere 
^4'^ioni respectively). It should be noted that the table entries corresponding 
to non-dispersive errors are zero, as a result of the multi-frequency calibration. Table [2] in¬ 
cludes also the dynamic range (DR) values measured from the SFPR maps, which are used 
to estimate a per baseline thermal phase error {(Jfpihm) using the expression 360°/DR/\/Nant, 
where Nant is the number of antennas. This is derived using the relation ship between po¬ 
sitiona l error and dynamic range {ao,,s ~ ^beam/DR) and the formulae in [Thompson et ah 


(120071. A12.58). Note that the dominant error contribution in Table l2]is related to the static 
component of the ionospheric propagation, which reaches peak values for frequency pairs 
with = 22 GHz and larger values of R and source pair angular separations; this will be 
revisited in the discussions section. The quadratic sum of the forementioned errors (\/Scr^) is 
converted to the final SFPR astrometric error (ctaocos^, cta^), for the KVN baselines lengths 
of ca. 400 km. 


Finally, we convert the SFPR astrometric errors to frequency dependent position shift 
errors for each source by: 1) passing those through the same SVD transformation used for 
the decomposition of the measurements, and 2) combine the outcome with the errors in the 
global shift minimization analysis, i.e. the errors in the jet position angles measured from 
the maps, as listed in Table [31 The hnal astrometric accuracies are listed in Table [31 in /uas, 

as O'/\a cos 6 Und 0/\^. 

These correspond to the errors in the measurements of the position shifts that enable the 
bona fide astrometric registration of the maps across frequencies, for each of the five observed 


sources. 
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5. Discussions 

5.1. Demonstration of Multi-Frequency Calibration and Astrometry up to 

130 GHz 


Simultaneous multi-frequency observations offer an effective path to achieve increased 
sensitivity and precision astrometry in mm-VLBI, beyond the domain of standard techniques, 
such as phase referencing. The compensation of the fast phase changes imposed by the rapid 
tropospheric fluctuations in mm-VLBI, using observations at a lower frequency of the same 
source, results in an increased coherence time of up to 20 minutes at 130 GHz using FPT 
analysis, which results in a significant increase of sensitivity. Moreover, when combined 
with the observations of another source, a bona fide astrometric measurement of the relative 
frequency dependent position-shift between the two frequencies can be estimated using the 
SFPR technique. This in turn results in an unlimited extension of the coherence time. 


The work presented in this paper corresponds to a hrst demonstration of SFPR at 
130 GHz, the highest frequency of the KVN. The application of SFPR techniques has al¬ 
lowed the detection of weak sources that were not directly detected within the atmospheric 
coherence time (i.e. with self calibration) and we have measured frequency dependent posi¬ 
tion shifts between a range of frequencies from 22 up to 130 GHz with high precision, for each 
of the observed AGNs. Previous attempts to carry out astrometry at such frequencies, with 
the VLBA up to 87 GHz, were very limited: Once with conventional phase r eferencing using 
very fast source switching and a very close source pair with ^IT separation flPorcas fc Rioia 


20021 1 and once with fast frequency switching using SFPR flRioia Sz DodsonI 1201111 . More¬ 


over, the fast frequency switching observing mode of the VLBA leads to residuals in the 
tropospheric compensation, which ultimately limit the accuracy and quality of the analysis. 
The simultaneous multi-frequency observing capability simplihes and widens the application 
to even higher frequencies and to many targets. Therefore we have demonstrated the ben¬ 
efits of simultaneous multi-frequency observations for mm-VLBI, to achieve sensitivity and 
astrometry at frequencies up to 130 GHz, the maximum frequency available with KVN. We 
believe these benefits would continue to apply beyond this frequency. 


5.2. Interpretation of the Measnred Position Shifts 

The outcome of the astrometric analysis presented in this paper is a measurement of 
the frequency dependent absolute position shift of the brightest features, or reference points, 
in the maps at the four KVN frequency bands, for each of the five AGNs. In general, the 
reference points in the KVN maps do not correspond to the position of the “core” compo- 
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nent, due to structure blending effects resulting from the relatively short ~400 km baselines. 
Therefore, in general, our measurements correspond to the position shifts or angular sepa¬ 
rations between the centroids of the brightness distributions at the four frequencies; in the 
case of a point source, this would be the same as the “core-shift”. 


Higher resolution observations with longer baselines would make it possible to isolate 
the “core” component as the reference point in the astrometric analysis, and to achieve 
an increase in the astrometric precision directly proportional to the enlarged baseline. For 
example, a 8000 km baseline would result in a twentyfold decrease of the astrometric errors 
listed in Tables [2] and |5l which would be suitable to measure the small magnitude of the 
opacity “core-shift” effect between the four KVN freque ncy bands predicted in the standard 
model for extragalactic radio sources flJung et all 120151) . Regardless of the baseline length, 
the SFPR measurements provide a bona hde astrometric registration of the maps, which 
is at the base for reliable spectral index maps and in general spectral distribution studies. 
Applying the absolute single source postition shifts in Table |5] to the hybrid maps in Figure 
|3] provides the required astrometric image registration for such analysis. 


In this paper, despite having the measurements of the shifts for the “bona hde” astro- 
metric registration of the maps, the poor amplitude calibration in our observations prevented 
us from obtaining meaningful spectral index maps. In a second epoch of observations we 
have included an improved amplitude calibration strategy to overcome this issue. 


The canonical SFPR astrometric errors range from a few tens to a few hundreds of 
pas depending on the frequency pairs, and are completely dominated by the systematic 
static ionospheric terms (erosion)- If this was actually the case we would expect to hnd 
(i) very similar SFPR astrometric offsets in all source pairs with similar separations, for a 
given frequency pair (i.e. 1928-1-738 and 2007-1-777 to 1803-1-784) and (ii) predictable ratios 
between the astrometric offsets, e.g. the ratio of the astrometric offsets for K—)-W and Q—?-W 
should be equal to 5 if the ionospheric uncertainties were dominant. As we see no indication 
of these signatures in our results we believe that the ionospheric contributions are acting 
coherently across the sources and that the canonical errors are therefore overestimates. This 
is quite possible as the KVN baseline lengths are of the same order as the height of the 
ionospheric E-layer, therefore the atmospheres over the different antenna sites would not 
be fully decorrelated. For this reason, we used an uncertainty in the total electron content 
(TEC) parameter in our error analysis equal to 3 TEC units (TECU), which we believe is 
still an over-estimate. 


In our results we have presented an experimental demonstration of successful “bona 
hde” astrometric registration between mm-VLBI maps, at the four KVN frequency bands, 
including, for the hrst time, high precision astrometric measurements at 130 GHz. The 
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combi nation with similar studies at lower frequencies fsee iRos et al.l (120011) : iGuirado et al. 


(120001) for wide-held high precision astrometry of some of these sources) opens a promising 
path towards providing a complete picture of the underlying physical mechanisms of jet 
formation across a very wide frequency range. 


5.3. Comparison of SFPR to other Methods 


There exist a variety of methods used to register the images of AGNs and/or maser 
species at different f r equen cies; here we discuss those, in comparison to SFPR, for mm- 
VLBI. iHovatta et ahl (120141) provides a thorough review of the alternative methods for reg¬ 
istration of the images of AGNs at m ultiple frequencies. These broadly consists of: using 
optically thin bright jet components (IFromm et al.l 120131) . 2D cross-correlation algorithms 
(IGroke fc Gabuzdal 120081 ). or a combination of both. However these methods are predicated 
on the assumption that there is a clearly identihable optically thin bright jet component, 
which can act as a reference point for all frequencies, or an ensemble of less bright optically 
thin jet components, which can provide an average registration. Hovatta et al. carried out 
an error analysis of the propagation of incorrect alignments on spectral index and rotation 
measure maps. They concluded that the alignment errors are dominant around the core 
region (up to a distance of ~3 mas) and therefore the conclusions on the spectral index dis¬ 
tribution for the innermost jet regions should be treated with caution. The application of 
these methods for compact sources and for sources with faint or smooth jets is clearly an 
issue. Therefore these methods can be unreliable or impossible to use in mm-VLBI, where 
in many cases only the compact core can be detected. 

In a similar fashion, some maser species can be assumed to form in a ring, and the centre 
of the ring can therefore act as a reference point (jPesmurs et al.ll2000[) . or one can identify a 
single component which appears similar in velocity and orientation with respect to the mai n 
body of emission and use that as the reference point across frequencies (iGotton et al.l 120041) . 
It is not hard to see the short-falls in such approaches and these different methods tend to 
produce incompatible conclusions. Phase referencing would provide a clear solution for such 
challenges but, as pointed out previously, phase-referencing is not an option above 43 GHz, 
in general. 


Therefore SFPR stands alone as a method which will allow for the unambiguous regis¬ 
tration across wide frequency spans for mm-VLBI images, both for continuum and spectral 
line studies, because it provides a complete compensation of atmospheric propagation and 
instrumental effects. SFPR is widely applicable for many sources, since the calibrator source 
can be at a signihcant angular separation and slow source switching does not undermine the 






















result. The method will work even at very high frequencies, making it particularly suitable 
for mm-VLBI observations. We do note that systematic effects do need to be carefully taken 
into account, particularly when using a lower frequency of 22 GHz. 


6. Conclusions 

We have demonstrated that the KVN multi-band system is capable of delivering in¬ 
creased coherence times by calibrating the highest frequencies with the scaled up phase so¬ 
lutions from the lower frequencies, using the Frequency Phase Transfer method. At 130 GHz 
the coherence times were extended from a few tens of seconds to 20 minutes, and to many 
hours by interleaving observations of a second source. This provides improved sensitivity 
through allowing longer integrations on weak sources. 

We have demonstrated that the KVN multi-band system is capable of delivering as¬ 
trometric results at the highest frequencies, using the Source Frequency Phase Referencing 
technique. We have measured accurate relative position shifts between frequencies in the 
range of 22 to 130 GHz, using observations of six pairs of sources with angular separations 
between 3.6° and 11°. 

We have shown how to decompose these relative measurements into absolute single 
source frequency dependent position-shifts. These decomposed position-shift measurements 
are all that is required to form high fidelity spectral index maps between the four frequency 
bands, which we will present in a subsequent paper. 
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